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Abstract 

> 

(f-} ■ In this work we essentially reinterpreted the Sieczka-Holyst (SH) 

in : 

fN) . model to make it more suited for description of real markets. For 

\ instance, this reinterpretation made it possible to consider agents as 

m 

crafty. These agents encourage their neighbors to buy some stocks if 



X 



agents have an opportunity to sell these stocks. Also, agents encour- 



H ■ age them to sell some stocks if agents have an opposite opportunity. 

Furthermore, in our interpretation price changes respond only to the 
agents' opinions change. This kind of respond protects the stock mar- 
ket dynamics against the paradox (present in the SH model), where 
all agents e.g. buy stocks while the corresponding prices remain un- 
changed. In this work we found circumstances, where distributions of 
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returns (obtained for quite different time scales) either obey power- 
law or have at least fat tails. We obtained these distributions from 
numerical simulations performed in the frame of our approach. 

PACS numbers: 89.65. Gh, 89.75.Fb, 02.50.Ey 



1 Introduction 



The market modeling is a modern challenge, especially in the context of still 
lasting a global financial crisis. Besides, the market modeling is a canoni- 
cal branch of econophysics lH3|. One of the significant part of this branch 
is agent based modeling inspired by Ising and Potts models. Notably, the 
approach based on the Potts model and its different generalizations are par- 



ticularly fruitful 4HlO|. On this way already several, although not all styl- 
ized facts, are well reproduced, e.g. fat tails of returns' distribution as well 
as long-range correlations of absolute returns, and short-range correlated re- 
turns. Certainly, this is a significant encourage to deal with these models. 
Nevertheless, since two decades of fast development of econophysics, there 
is still no definite model reproducing all known up to now stylized facts 
concerning real financial markets. Furthermore, even a microeconomic (or 
microscopic) sources of some stylized facts are still not well understood. 

In this paper we mainly reinterpret the Sieczka-Holyst (SH) threshold 
model of financial markets Jllj], taking into account the concept of the ne- 
gotiation round from the Iori model 6|. We choose the SH model for the 



following reasons: 

(i) It contains a threshold mechanism with three possible values of spin 
variable. 

(ii) It considers interaction between agents as well as noise in their activity. 

(iii) It assumes strength of the interaction which strongly depends on the 
macroscopic state of the system. 

(iv) The SH model seems to be sensitive to the concrete kind of emotions 
which each agent subordinates during the stock market evolution 
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(v) The SH model seems to be a more realistic than other agent-based 
models. 

Despite of these properties, the SH model leads to the paradoxal situation 
where if all agents buy stocks or all sell them, the corresponding prices remain 
still unchanged. This paradox inspired us for the appropriate reinterpretation 
of the SH model. 

Moreover, in our approach we model emotions of agents, including even 



very intensive ones, by the mean of the Weierstrass-Mandelbrot noise 
instead of the Gaussian noise considered in the SH model. 
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2 Sketch of the SH model 



Let us shortly remind the SH model. There are N interacting agents placed 
on 2-dimensional square lattice nxn, where N = n 2 . Agent number i (where 
2 = 1,..., N) can take into account only one of the three actions: buy, stay 
inactive or sell. Such an agent is represented by 3-state spin variable Sj with 
possible values: +1 if agent is buying stocks, if he stays inactive, and — 1 
in the case of selling stocks. 

The single time step t, called a round, consists of N drawings of spinal . 
Hence, in a given round each agent has, in average, a single chance to change 
his opinion. After each drawing the chosen spin is updated according to the 
rule 

Si (t) = sign A | M(i _ 1} 



N 



3=1 



1) 



where 



+ 1 if x > q, 

if -q < x < q, (2) 
— 1 if x < —q, 

and constants A, a are positive. The value of pair interaction strength = J 
if agent j is one of the four nearest neighbors of agent i, otherwise J^- = 
0. The temporal noise r)i(t) is a random variable drawn from the standard 
normal distribution. It represents individual erratic opinion of i-th agent. 
Furthermore, the change of any spin can affect its neighbors immediately 



1 Single round resembles lMCS/spin used in dynamic Monte Carlo methods. 
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(i.e. within the same round). 

The magnetization of the network, present in Expression ([1]) is defined as 
usual, 

1 N 

M(t) = -Y,s l (t). (3) 

The crucial role of the magnetization, as a quantity determining activity of 
agents, is well seen. That is, index X\M(t — 1)| in the sign function is a 
temporal threshold. Hence, if the absolute value of the neighbors' impact on 
the i-th agent together with his (erratic) opinion is lower than the threshold, 
the agent does not participate in the market game assuming a neutral position 
(as his spin variable equals 0). In other words, the large absolute value of 
magnetization gives (at a fixed value of A) the large value of the threshold. 
Hence, the probability that each agent assumes a neutral position distinctly 
increases. 

As the value of magnetization was already defined, the temporal price of 
a stock 

P{t) ocexp(M(*)). (4) 

In the simulation Sieczka and Holyst used a square lattice of agents with 
n = 32 and N = 1024, possessing the periodic boundary conditions, and 
assuming the initial configuration of spins as random. For thermalization of 
the system they used the first 5000 rounds. 

As authors claimed, simulations' results reproduced main stylized facts 
such as volatility clustering in the variograms of returns, fat (but not the de- 



5 



sired power-law) tails of returns' statistics and fast decaying autocorrelation 
function of returns revealing a negative feedback. However, autocorrelation 
of absolute returns decayed exponentially and not according to the desired 
power-law. In the frame of our approach we repaired both the statistics of 
returns and autocorrelation of absolute returns. 

3 Reinterpretation of the SH model 

In this Section, we reinterpret the SH model by essential change of the mean- 
ing of the spin variable Sj (i = 1, . . . , N). We postulate that instead of an 
action of the agent, this variable means only the agent's opinion about the 
market situation, which he communicates to his nearest neighbors in ac- 
cordance with Expression ([I]). That is, Sj = +1 means a positive opinion 
communicated by the agent (i.e. according to this message his nearest neigh- 
bors should buy stocks). For = -1 a negative opinion is communicated 
by him (i.e. his nearest neighbors should sell stocks). Value Sj = simply 
means a lack of the agent opinion. 

In our approach the activity of the agent always requires two subsequent 
opinion stages. That is, the change of spin variable di(t) = Si(t) — Si(t — 1) 
(i.e. the change of the agent's opinion during the subsequent rounds) means, 
for di > 0, the agent's demand and reversely, for d\ < 0, this change means 
the agent's supply (considered here as a negative demand). 

In other words, the agent declares demand for stocks if his spin variable 
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increased within a given round (in comparison with value of the spin variable 
obtained within the previous round), while he offers a supply of stocks if the 
value of his spin variable decreased. For instance, the change of the spin 
variable from —1 to +1 (or from +1 to —1) defines the largest possible 
demand (or supply) declared by a given agent in a single round. This is 
because the absolute value of the difference between subsequent spin values 
is the largest one, i.e. = 2. Other changes define the smaller demand 

or supply as then = 1. We assume that declared demand or supply 

of any agent is realized immediately, i.e. in the same round. As a result, 
the agent buys or sells stocks. Note that in this approach the agent opinion 
concerning the market situation always remains unchanged (by definition) 
until the next transaction will occur. Thus, two opinion stages (i.e. two 
subsequent different values of agent's spin) were realized. 

Obviously, after buying some stocks the agent immediately (by definition; 
cf. Expression (TTJ) communicates the optimistic or neutral opinion to his 
neighbors, as he is interested in an increase of a stock price. In the opposite 
case (i.e. in the case of selling), the agent gives a pessimistic or neutral 
opinion to his neighbors (as he is interested in a decrease of a stock price 
directly after the selling). This is the way how the subsequent two-stage 
cycle of the agent activity (i.e. buying or selling and communication) is 
repeating. 

One of the consequences of our reinterpretation of spin variable is that 
the agents are assumed crafty. That is, if the spin value of agent = +1, he 
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communicates positive opinion to his neighbors but he cannot buy a stock 
by himself. It is because his spin value cannot increase, i.e. he can only keep 
or sell stocks. For Sj = — 1 the situation is reversed, i.e. the agent which 
earlier sold some stocks, now can only keep or buy them. In this way, the 
negative coupling between subsequent single agent's opinions can be realized. 
For Si = the agent occupies a neutral position and he can either buy or sell 
a stock as well. Indeed, in this paragraph the principal differences between 
our approach and the SH model were precisely defined. 

In our approach the noise term crr)i(t) indicates an intrinsic (erratic) opi- 
nion of agent i and parameter a is a strength of this opinion. It should not 
be mixed up with the spin variable Si(t), which represents opinion directly 
communicated to other agents. 

Our reinterpretation of the spin variable protects the stock market dy- 
namics against the situation (present in the SH model), where all agents buy 
stocks while corresponding temporal prices remain unchanged (as magneti- 
zation becomes unchanged). Furthermore, our approach ensures correlation 
between volatility of returns and a trading volume, as expected. 

As another significant consequence of our approach, the temporal thresh- 
old is able to slow down both the fast increase of the price if trend is positive 
and its fast decrease if trend is negative. That is, the threshold introduces 
a kind of a price damping, which cannot lead neither to oscillations nor to 
change of trend direction. 

Usually, in agent-based models one considers simplified linear model of a 
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price formation, where change of price (or logarithmic return) is proportional 
to excess demand [15] . Namely, 



lnP(t) - In P(t-r) 



dcf. 



r T (t) oc ED T {t). 



(5) 



Hence, for the particular case r = 1 the excess demand can be written in the 
following form, 



Although Expression ([5]) is formally equivalent to Equation fl4]), it has a 
different interpretation as spin Si(t) means here an opinion of i-th agent and 
not his individual demand (as it is in the SH model). Apparently, the excess 
demand, ED(t), can change if and only if the opinion of any agent changes. 
Notably, in the SH model the excess demand does not change although an 
agent still buys a stock. Due to our reinterpretation of a spin variable, we 
avoided this paradox. 

In the SH model the noise is significant only to some extent. That is, if 
the value of the noise exceeds the range [—(4 J + A), 4 J + A], the exact value 
of this noise is then unimportant. This comes from the threshold character 
of Expression ([T|) . 

In Expression ([T]), there are present three parameters which (without loss 
of generality) can be reduced to two crucial relative parameters: 



N N 



ED(t) = £ di(t) = 5>«(*) -Si(t- 1)). 



(6) 



i=l i=l 
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(i) parameter a.\ — jj. If ax < 1, then the nearest neighbors' impact on 
the agent is high, otherwise the agent is an opportunist or skeptical as 
he often chooses Sj = 0. 

(ii) Parameter a a = -fj. If a a < 1, the agent trusts more his neighbors 
than himself, otherwise he trusts more himself than his neighbors. 



In paper ll| it was considered, for instance, the case ot\ > 1 and a a < 1, i.e. 
opportunist or skeptical agents, trusting more their neighbors than them- 
selves. If agents are opportunist or skeptical, then even coherent feedback 
coming from all neighbors (having the same value of spin variable) can leave 
the agent's opinion unchanged. It seems to be useful to systematically ex- 
amine all ranges of parameters ot\ and a a . 



4 Results and discussion 

We begin our simulation with randomly oriented spins distributed over the 
square lattice consisting of N = 1024 sites. In our simulation we used noise 
distribution in the form of the Weierstrass-Mandelbrot probability density 
function, 

/ 1 \ 00 1 1 
P(x) = (l - E xJ ■ " W> A '' b>l,b o >0. (7) 

We used this distribution for the following reasons: 

(i) Agent's emotions, which can contain even very intensive ones, should 
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be ruled by power-law distributions which may be better suited than 
canonical Gaussian ones. 

(ii) Herein, the control of the noise by its variance and exponent is possible. 
This variance can vary in the range from finite to infinite values. 

The Weierstrass-Mandelbrot distribution is a canonic simple one poss- 
esing above given properties. This distribution is discrete having spikes at 
x = ±Oo6 J , 3 — 0, 1, 2, . . .. Its variance is finite only for b 2 /K < 1 and takes 
the form 

9 ,, 1 - l/K 
<j = bl ^ 



°1 -b 2 /K v ' 

Otherwise, it is an infinite quantity. 

It can be easily proved that for |x| ^> 1/lno Expression (J2J) can be ap- 
proximated by 

. , I — l/K Oq , n 111 K 

pix) = — - — — — i , where exponent p = - — — . (9) 
lniv mo 

Note that a 2 is finite only for /3 > 2 otherwise, it is infinite. 

Although Equation (jHJ) is a power-law distribution, only the part of this 
distribution extended over the maximal range [— (4J+A), 4J+A] is significant 
for the system dynamics (as it is suggested by Expression ([TJ). Inside this 
maximal range spin can assume values equal to —1,0,-1-1. Above the right 
border of this range only value of spin equal to +1 appears while below the 
left border only —1 appears. How this affects the the power- law tail of the 



11 



noise distribution is still an open question. 

In Fig. [1] we plotted returns' histograms obtained in our simulation (in 
lin-log scale) for different values of time lag r. In the simulation we used 




Figure 1: Histograms of returns r T (t) from our simulation, for instance, for 
parameters J = 1,A = 1, K = 5, b = 2, bo = 0.21, and different time lag 
t. The returns were re-scaled by the corresponding standard deviations of 
returns' time series. The distributions were separated for the better view. 
The dashed line is a normal distribution fitted to the central part of histogram 
for r = 256. The solid line denotes the power-law with exponent equals 
3.322 = 1 + (see text for details). 



assumptions considered in Section [2j Besides, we set, for instance, J = 
1,A = 1, K — 5, b — 2, bo — 0.21. Hence and from Formula (jHJ), we obtain 
a = 0.42 and a a = 0.105 while ot\ = 0.25. We verified that for assumed 
parameters only components with j = 0, 1, 2, 3, 4 are significant values in the 
sum in Expression (j7j). Apparently, the histograms for the most (i.e. for 
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small and intermediate) values of r have well-formed power-law tails with 
exponent which (to good approximation) equals 1+/3, where — In Kj In b = 
In 5/ In 2 = 2.322. That is, these histograms quite well reproduce statistics 
given by Equation ([9]). Remarkably, although the noise is cut off by borders 
— (4 J + A) and 4 J + A, the histograms of returns have still power-law tails 
driven by exponent 1 + 0. The more so, a,\ < 1, which means that the range 
[—(4 J + A), 4 J + A] is narrower than the corresponding one in the SH model. 
Nevertheless, for the largest values of r the returns' distributions still exhibit 
fat tails. These results distinguish our approach from the SH model. 

In addition we explain herein how we avoid a possible trapping in our 
simulation. By trapping we understand such a single realization of the simu- 
lation, where after sufficiently long time all spins are spontaneously ordered 
(a fully ferromagnetic state). Then the market is trapped by a stable state 
of extreme magnetization (except the case of a\ > 1 or a a > 1 where this 
trapping does not occur at all). To avoid this trapping the system was acti- 
vated by some exogenous factor making an abrupt transition of the system 
to the paramagnetic statj^ (where all spins are randomly oriented). Then, 
the analysis of the system between these subsequent transitions is made. 
2 Abrupt transitions are a characteristic feature of modern financial markets. 
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5 Conclusions 



In this work we essentially reinterpreted the Sieczka-Holyst threshold model 
of financial markets. This reinterpretation means that we introduced a new 
meaning of the spin variable, although the mathematics of the SH model re- 
mained unchanged. In our approach the value of spin variable does not mean 
an action of the agent but, instead, it is only the agent's opinion about a mar- 
ket situation. This opinion he communicates to his neighbors. The agent's 
action (i.e. buying, selling or staying inactive) is in our approach connected 
only with the change of spin variable. Indeed, this protects the stock market 
dynamics against the paradox (present in the SH model), where all agents 
buy stocks while corresponding temporal prices remain still unchanged. 

Furthermore, our approach differs from the SH model as we used the 
Weierstrass-Mandelbrot noise to describe an intrinsic (erratic or even very 
intensive emotional) opinions of agents instead of the Gaussian noise used in 
the SH model. 

Despite using of a threshold, the results obtained for Gaussian and 
Weierstrass-Mandelbrot noises are quite different. The difference between 
both noises within the range [—(4 J + A), 4 J + A] has a significant impact 
on the results. For instance, for the Weierstrass-Mandelbrot distribution 
there are power-law statistics for (almost) full range of returns, while for the 
Gaussian one the corresponding range is too short to say something definite. 
Obviously, the discreteness of Weierstrass-Mandelbrot distribution influences 
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the results. 

We simulated the above sketched dynamics of agents and calculated 
statistics (histograms) of returns. Some of these statistics have power-law 
tails while others have only fat tails. Explanation of this result is still a chal- 
lenge. This is because power-law tails of the noise is cut-off by borders equal 
to ±(4J + A). Nevertheless, power-law tails of returns' statistics is driven by 
the same exponent /3 which governs the noise distribution. The systematic 
comparison of our approach predictions (additionally concerning e.g. auto- 
correlation functions and power spectrum) with empirical data would be a 

n 

decisive verification of this approach_|. 
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